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Abstract. We present fully self consistent 3D simulations of compressible Hall MHD plasma that describe spectral features 
relevant to the solar wind plasma. We find that a kT 1 ^ spectrum sets in for the fluctuations that are smaller than ion gyro 
radius. We further investigate scale dependent anisotropy led by nonlinear processes relevant to the solar wind plasma. Our 
work is important particularly in understanding the role of wave and nonlinear cascades in the evolution of the solar wind, 
structure formation at the largest scales. 
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1. INTRODUCTION 

Solar wind (SW) fluctuations comprise a multitude of 
length and time scales that collectively exhibit numer- 
ous nonlinear processes. Despite it's complex evolution- 
ary dynamics, solar wind provides the best laboratory 
for testing many nonlinear theories and simulation mod- 
els. Spacecraft databases provide compelling observa- 
tional evidences that the inertial range magnetic field 
fluctuations associated with characteristic frequencies, 
that are smaller than the ion gyro frequency (i2 n ), can 
be described predominantly by a Kolmogorov-like 5/3 
spectrum Theoretic analysis describing the 5/3 SW 
spectrum regime relies largely on the usual magnetohy- 
drodynamic (MHD) model of plasma. Intriguingly, the 
higher time resolution databases of solar wind fluctu- 
ations depict a spectral break near the end of the 5/3 
spectrum that corresponds to a high frequency (> Q ci ) 
regime where turbulent cascades are not explainable by 
the usual MHD models. This refers to a second inertial 
range where turbulent cascades follow a £~ 7 / 3 [2] (where 
k is a wavenumber) spectrum in which the characteris- 
tic fluctuations evolve typically on kinetic Alfven time 
scales. The onset of the second or the kinetic Alfven in- 
ertial range is still elusive to our understanding of SW 
turbulence and the issue has been under a constant de- 
bate since many years. The mechanism leading to the 
spectral break has been thought to be either mediated by 
the kinetic Alfven waves (KAWs) Jj], or damping of ion 
cyclotron waves, or dispersive processes [3], Hall effects 

In this paper, we describe results from our three di- 
mensional simulations that explains that the fc~ 7 / 3 spec- 
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trum observed above the spectral break in SW may be 
led by the Hall effects in the KAW regime. The un- 
derlying model, in our simulation model, is based on 
a two fluid Hall MHD description of plasma that con- 
sists of both electrons and ions. In the high frequency 
regime, co > £l c i, the inertialess electrons contribute to 
the electric field which is dominated essentially by the 
Hall term corresponding to J x B force. The latter, upon 
substituting in the ion momentum equation, modifies ion 
momentum, magnetic field and total energy in a man- 
ner to introduce a high frequency (co > Cl C i) and small 
scale (k±pi > 1, where p^ is ion Larmour radius) plasma 
motion. The characteristic length scales (k~ l ) associ- 
ated with the plasma motions are smaller than ion gyro 
radii (pi)- The quasi-neutral solar wind plasma density 
(p), velocity (U), magnetic field (B) and total pressure 
(P = P e + Pj) fluctuations can then be cast into a set of 
Hall MHD equations as follows (see Ref Ja] for detail 
description). 

^ + V-(pU) = 0, (1) 

p ^+U-V^U=-VP+^JxB + ju (Vu+^VV-U 

(2) ' 

^ = Vx (\]xB-di^^j +TjV 2 B, (3) 

de (\ , yP c \ 

3-+V- l-pU 2 U + -!—U + — ExB =0 (4) 

dt \2 y-l 4x ) 

where e = \/2pU 2 + P/(y- l)+B 2 /8n is a total en- 
ergy of plasma that contains both electron and ion mo- 
tions. All the dynamical variables are functions of three 
space and a time, i.e. (x,y,z,t), co-ordinates. Equa- 
tions (03 to (0]i are normalized by typical length £q 
and time fo = W v 'o scales in our simulations, where 
v'o = B()/(Anp(f) l l 2 is Alfven velocity such that V = 
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FIGURE 1. Inertial range turbulent spectra for magnetic and 
velocity field fluctuations. The fluctuations closely follow re- 
spectively fc~ 7//3 scalings in the kd, > 1 KAW regime. 

£ V, d/dt = t d/dt,lJ = U/v ,B = B/vo(47rp ) 1 / 2 , J P = 
P/PqVq,p = p I po. The parameters [I and 77 correspond 
respectively to ion-electron viscous drag term and mag- 
netic field diffusivity. While the viscous drag modifies 
the dissipation in plasma momentum in a nonlinear man- 
ner, the magnetic diffusion damps the small scale mag- 
netic field fluctuations linearly. The magnetic field is 
measured in the unit of Alfven velocity. The dimension- 
less parameter in magnetic field Eq. (O i.e. ion skin 
depth dj = dj/lo,dj — C/co P i is associated with the Hall 
term. This means the ion inertial scale length (di) is a 
natural or an intrinsic length scale present in the Hall 
MHD model which accounts for finite Larmour radius 
effects corresponding to high frequency oscillations in 
kdi > 1 regime. Clearly Hall forces dominate the mag- 
netoplasma dynamics when l/p(JxB)>UxB term in 
Eq. ([3]l which in turn introduces time scales correspond- 
ing to the high frequency plasma fluctuations in kdi > 1 
regime. Furthermore, our model includes a full energy 
equation [Eq. ©J unlike an adiabatic relation between 
the pressure and density. The use of energy equation en- 
ables us to study a self-consistent evolution of turbulent 
heating resulting from nonlinear energy cascades in the 
solar wind plasma. 

2. TURBULENT SPECTRA 

To study the nonlinear evolution of turbulent cascades 
in a Hall MHD solar wind plasma, we have developed 
a fully three dimensional compressible Hall MHD code. 
Our code is massively parallelized using Message Pass- 
ing Interface (MPI) to run on cluster like distributed- 
memory supercomputers. The code is scalable and trans- 
portable on different cluster machines. The spatial des- 
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FIGURE 2. Evolution of (k±(t)) and (kn(t)) as a measure 
of anisotropy in Hall MHD. Initially, (k ± (t = 0)) = k» (t = 0). 
As time evolves, anisotropy in spectral transfer progressively 
develops such that (k±(t)) > (ku(t)}. 

critization employs a pseudospectral algorithm [7] based 
on a Fourier harmonic expansion of the bases for phys- 
ical variables (i.e. the density, magnetic field, velocity, 
temperature and energy) whereas the temporal integra- 
tion uses a Runge Kutta (RK) 4th order method. The 
boundary conditions are periodic along the x,y and z di- 
rections in the local rectangular region of the solar wind 
plasma. 

The turbulent fluctuations are initialized by using 
a uniform isotropic random spectral distribution of 
Fourier modes concentrated in a smaller band of lower 
wavenumbers (k < 0.1 k max ). While spectral amplitude of 
the fluctuations is random for each Fourier coefficient, it 
follows a certain initial spectral distribution proportional 
to k~ a , where a is initial spectral index. The spectral dis- 
tribution set up in this manner initializes random scale 
turbulent fluctuations. A constant background magnetic 
field is included along the z direction to deal primarily 
with the large scale or background solar wind magnetic 
field. Turbulent fluctuations in our simulations are driven 
either at the lowest Fourier modes or evolve freely under 
the influence of self-consistent dynamics described by 
the set of Eqs. (Q]) to (0). The inertial range spectral cas- 
cades in the either cases lead to the nearly identical turbu- 
lent spectra. We have further carried out simulations for 
a range of various parameters and spectral distributions 
to ensure the validity of our codes and the physical re- 
sults. The simulation parameters are; spectral resolution 
is 128 3 , 77 = (X = 10 3 ,j8 = 1.0, kdi - 0.1 - 10, L x =Ly = 
L z = 2n. The nonlinear coupling of velocity and mag- 
netic field fluctuations, amidst density perturbations, ex- 
cites high-frequency and short wavelength (by the 0)/(B 6 , 
effect) compressional dispersive KAWs. The nonlinear 
spectral cascade in the modified KAW regime leads to 
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FIGURE 3. Spectrum of and fey as a function of k. 
The large scale inertial range turbulent fluctuations are more 
anisotropic as compared to the smaller ones. 



a secondary inertial range in the vicinity of kdi ~ 1, 
where the turbulent magnetic and velocity fluctuations 
form spectra close to kr 1 ^ Jfl [gl [l(J . This is displayed 
in Fig. ©. It is shown in Shaikh & Zank (2009) that the 
spectra described in Fig (1) is led by Hall effects. 

3. ANISOTROPIC CASCADES 

To measure the degree of anisotropic cascades (or 
spectral anisotropy) by employ the following diag- 
nostics to monitor the evolution of the k± mode in 
time. The averaged k± mode is determined by av- 
eraging over the entire turbulent spectrum weighted 

by £ ± ,thus (*j.(t)) = (Lt|*j.C2(*,t)l 2 /EtlC2(*,0l 2 ) 1/2 - 

Here (•••) represents an average over the entire Fourier 
spectrum, k^ = (k 2 + k 2 ) l l 2 and Q represents any of 
B,V,p, V x B and V x V. Similarly, the evolution of 
the k\\ mode is determined by the relation, (&||(f)) = 

(Lk\k\\Q(k,t)\ 2 /Lk\Q(k,t)\ 2 y/ 2 . It is clear from these 
expressions that the (k±(t)} and modes exhibit 

isotropy when (k±(t)) ~ (kn(t)). Any deviation from 
this equality corresponds to spectral anisotropy. We fol- 
low the evolution of the (k±(t)) and (kn(t)) modes 
in our simulations. Our simulation results describing 
the evolution of (k±) and (k\t(t)) modes are shown in 
Fig. It is evident from Fig. (f2]i that the initially 
isotropic modes (k±(t)) ~ (kn(t)) gradually evolve to- 
wards a highly anisotropic state in that spectral transfer 
preferentially occurs in the (k±(t)) mode, and is sup- 
pressed in (fc||(f)) mode. Consequently, spectral trans- 
fer in the (k±(t)) mode dominates the nonlinear evolu- 
tion of fluctuations in Hall MHD, and mode structures 
become elongated along the mean magnetic field or z- 



direction. Hence nonlinear interactions led by the nonlin- 
ear terms in the presence of background gradients lead to 
anisotropic turbulent cascades in the inertial range turbu- 
lent spectra. 

Fig. (f2]i illustrates anisotropy corresponding to an av- 
eraged k mode but the anisotropy exhibited by the small 
and large scale B and v fluctuations is not distinctively 
clear, nor is the degree of anisotropy in B and v fields 
clear from Fig. (0. The scale dependece of turbulent 
anisotropy is described in Fig. 0. Fig. (0 shows dis- 
crepancy in k± and k\\ is prominent at the smaller k's. 
This essentially means that the large scale turbulent fluc- 
tuations are more anisotropic than the smaller ones in 
a regime where characteristic length scales are smaller 
than dj i.e. kdi > 1- It further appears from Fig. Q that 
the smaller scales in the kdi > 1 are virtually unaffected 
by anisotropic kinetic Alfven waves that propagate along 
the externally imposed mean magnetic field Bq. Turbu- 
lent fluctuations with small characteristic scales in the 
kd, > 1 regime of Hall MHD are not affected by the 
mean magnetic field or kinetic Alfven waves. This leads 
us to conjecture that small scale turbulence in the kdi > 1 
regime behaves essentially hydrodynamically i.e. as ed- 
dies independent of the mean magnetic field or collision- 
less magnetized waves. Thus large and smaller turbulent 
length scales evolve differently in the kdi > 1 regime of 
Hall MHD turbulence. 



3. DYNAMICAL ALIGNMENT OF 
FLUCTUATIONS 

To understand the strength of the nonlinear interactions 
in Hall MHD solar wind plasma, we determine the de- 
gree of alignment of the velocity and magnetic field fluc- 
tuations by defining the following alignment parameter 
111 ill that spans the entire ^-spectrum in both the kdi > 1 
(Hall MHD) and kdi < 1 (usual MHD) regimes, ©(f) = 
cos- 1 (I k V k (0 -B k (0/Ek |V k (0||B k (0l) • The summa- 
tion is determined from the modes by summing over the 
entire spectrum. In this sense, the alignment parameter 
depicts an average evolution of the alignment of veloc- 
ity relative to the magnetic field fluctuations. Note care- 
fully that this alignment can vary locally from smaller 
to larger scales, but the averaging (i.e. summing over the 
entire spectrum) rules out any such possibility in our sim- 
ulations. Nonetheless, defined as above enables us to 
quantitatively measure the average alignment of the mag- 
netic and velocity field fluctuations while the nonlinear 
interactions evolve in a turbulent solar wind plasma. 

The Alfvenic cascade regime of MHD turbulence 
kdi < 1 in the solar wind plasma possesses relatively 
large scales (kdi < 1) in which the velocity and magnetic 
field fluctuations are observed to be somewhat obliquely 




FIGURE 4. A progressive decrease in the angle of alignment from 90° indicates the eventual weakening of V x B nonlinear 
interactions. Evolution of the degree of alignment in the kdj > 1 regime of Hall MHD. Small scale fluctuations possess orthogonal 
velocity and magnetic fields. 



aligned. Hence our simulations show that the angle of 
alignment evolves towards < 90°, as depicted in Fig. 
©. Hence the strength of the nonlinear interactions cor- 
responding to the V x B nonlinearity is relatively weak. 
This result is to be contrasted with characteristic turbu- 
lent length scales in the kdi > 1 regime. The angle of 
alignment for the smaller scales corresponding to the 
kdj > 1 regime is shown in Fig. ©. Significant differ- 
ences are apparent in the angle of alignments associated 
with the large and small scales Fig. It appears from 
our simulations that the small scale fluctuations (kdi > 1) 
are nearly orthogonal as seen in Fig. ©. By contrast, the 
large scale fluctuations (kdj < 1) in Fig. © show a sig- 
nificant departure from the orthogonality. 

4. CONCLUSIONS 

In summary, we have investigated the nonlinear and tur- 
bulent behavior of a two fluid, compressible, three di- 
mensional Hall MHD model. In the presence of a large 
scale mean background magnetic field, small scale turbu- 
lent fluctuations exhibit anisotropic power spectra close 
to £~ 7//3 . We find that the long length scales in the kdj > 
1 KAW regime are more anisotropic compared to the 
shorter scales. Dynamical alignment and angular distri- 
bution of turbulence velocity and magnetic field fluctua- 
tions is found to play a critical role in determining the de- 
gree of nonlinear interactions. We find that characteristic 
turbulent flucutations in the kdj > 1 regime relax towards 
orthogonality, so that most of turbulent scale fluctuations 
have velocity and magnetic fields that are nearly orthog- 
onal , i.e. V _L B. For large scale fluctuations correspond- 
ing to the MHD regime, magnetic and velocity fields are 
not perfectly orthogonal, being instead on average nearly 



70° to each other. By contrast, small scale fluctuations in 
the kdj > 1 KAW regime exhibit nearly perfect orthog- 
onality in that the average magnetic and velocity fields 
make an angle of nearly 90° with respect to each other. 

The spectral properties of nonlinear Hall MHD are 
particularly relevant for understanding the observed so- 
lar wind and heliospheric turbulence. Hall MHD may 
also be useful for understanding multi-scale electromag- 
netic fluctuations and magnetic field reconnection in the 
Earth's magnetosphere and in laboratory plasmas. 
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